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ABSTRACT 


Experimentally determined values of the local skin-friction coeffi- 
cients of adiabatic flat-plate, turbulent, compressible boundary layers 
were approximated by an analytical expression which, in turn, was 
used in conjunction with the boundary-layer integral-momentum 
equation to calculate test-section boundary-layer thicknesses in the 
Jet Propulsion Laboratory 20-in. Supersonic Wind Tunnel. These cal- 
culations were compared with measured boundary-layer thicknesses 
for the Mach number range of 1.4 to 5. These comparisons were fair, 
but varied as the Reynolds number was varied. The use of a simpler, 
analytical boundary-layer skin-friction equation gave better compari- 
sons and was not as dependent upon Reynolds number as the some- 
what less simple analytical expression. A short-cut method (in contrast 
to the step-by-step integration from the throat to the test section) foi 
estimating the test-section boundary-layer momentum and displace- 
ment thicknesses is shown, and several examples are given. 


I. INTRODUCTION 


Many experimenters have measured the local skin- 
friction of adiabatic flat-plate, turbulent boundary layers 
in the ranges of 1.5 < M < 5.8 and 10' < R„ < 15 X 10' 
(see Refs. 1-7). Since very careful measurements of the 
Jet Propulsion Laboratory (JPL) 20-in. Supersonic Wind 
Tunnel (SWT) test-section boundary layer have been 
made for 1.4 < M < 5 (mainly on the centerline of the 
two-dimensional nozzle curved walls; see Ref. 8) it 
seemed appropriate to determine the degree of compati- 
bility of these two types of measurements when related 
by means of the boundary -layer integral momentum 
equation. 


The procedure followed was: (a) to establish a simple 
analytical expression for the local skin friction from a 
turbulent boundary layer over an experimental flat plate, 
(b) to use experimental boundary - layer displacement - 
thickness to boundary -layer momentum -thickness ratios, 
and then (c) to apply the boundary -layer integral mo- 
mentum equation using the designed Mach-number dis- 
tribution along the curved walls of the JPL 20-m. SWT. 
This was done for test-section Mach numbers of 1.4, 1.64, 
2, 3, 4, and 5 at maximum tunnel supply pressures, and 
for M — 2, 3, and 4 at considerably decreased supply 
pressures. 
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II. EXPERIMENTAL LOCAL SKIN-FRICTION OF THE TURBULENT 
BOUNDARY LAYER OVER AN ADIABATIC FLAT PLATE 


Figure 1 presents the experimental local skin-friction of 
the turbulent boundary layer over an adiabatic flat plate 
(Refs. 1-7) as a function of the momentum -thickness 
Reynolds number. The skin-friction data were normalized 
by being ratioed to the corresponding Schultz -Grunow 
data on the local skin friction of an incompressible, turbu- 
lent boundary layer over an adiabatic flat plate (Ref. 9) 
as shown in Fig. 2. 

r _ 0.0334 

^ (log 10 R 9 ) 1-838 (1) 

The experimental data shown in Fig. 1 were faired with 
straight lines of equal slopes, which in turn were cross- 
plotted in Fig. 3, giving good agreement with the 
analytical, Mach-number-dependent, local skin-friction 
equation (Ref. 10). 


C, 

-jrT- = (1 + 0.144 A/2) -0.578 (2) 

The data in Figs. 2 and 3 were then approximated by the 
Mach- and Reynolds-number-dependent expression [Eq. 
(3)], the result being shown in Fig. 4. 



Fig. 1. Measured local skin-friction coefficient of the 
turbulent boundary layer over an adiabatic 
flat plate; C,JC,. vs. R„ 


C, 

= (1 + 0.144 M 2 ) °- 57s X 


^ 2 X 10^ " "* 


(3) 
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III. CALCULATION OF THE BOUNDARY LAYER IN THE JPL 20-IN. SWT 


The distribution of the Mach number along a curved 
wall of the JPL 20-in. SWT two-dimensional nozzle is 
shown in Fig. 5 along with the location of the downstream 
end of the curved wall. The boundary-layer momentum 
thickness at the nozzle throat was obtained by the method 
used in Ref. 11. 


<9 0 (ft) = 5.06 X 10~ 3 


(v*) 1 ' 5 X (t*L*) 2/5 
(T t ) 1/10 


(4) 


The values of 0„ used are listed in Table 1 along with 
data from Ref. 8. 


The one relation left to determine before being able 
to use the boundary-layer integral-momentum equation 
[Eq. (5)] is the variation of the boundary-layer displace- 
ment to momentum thickness ( H ). 

do _ C f dU w f 2 — M 2 -f H 1 

dx 2 dX X |_M (1 + 0.2M 2 ) y (5) 


H 0 = 1.3 + 0.46 M 2 (6)* 

The ratio H/H 0 increases with decreasing Reynolds num- 
ber, and remains relatively constant for M > 2.5. The 
total excursion of H/H 0 throughout the Mach number 
range, for the maximum supply pressure case, is from 0.96 
to 1.02, a total variation of only 6%. In the integration of 
Eq. (5), the value H/H 0 was made a constant, depend- 
ent upon each nozzle Mach number and the supply pres- 
sure condition (high or low). These values of H/H 0 were 
chosen to be somewhat of an average from the throat to 
the test section, and are listed in Table 1. The integration 
of Eq. (5) is not particularly sensitive to the chosen ratio 
of H/H 0> a 3% variation giving only a 1% change in the 
test section boundary-layer momentum thickness ( in 
the opposite direction) for the M = 3 nozzle. Of course, 
the boundary-layer displacement thickness will then be 
changed by 2% (in the same direction as the change in 
the H/H 0 ratio). 


The variation of test-section H with Mach number and, 
to a limited degree, with Reynolds number was obtained 
from the data of Table B-l in Appendix B (see Ref. 8; the 
value of H at the throat is included from Ref. 11). This 
variation, shown in Fig. 6, was normalized by ratioing it 
to the Hy value of Ref. 10, which is plotted in Fig. 7. 


The results of these wind-tunnel nozzle, boundary-layer 
thickness calculations are shown in Table 2, where they 
are compared with the measured values. As can be seen, 
the calculated values of the boundary-layer thicknesses 

a In air for y = 1.4. 
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Table 1 . JPL 20-in. SWT test-section boundary-layer measurements 
on centerline of curved wall 


M 

Pt 

h 

R/ in. 

X 

t * 

c 

A 

e m 

Hm 

c b 

H* c 

cm Hg 

°R 

X 10' 5 

in. 

in. 

in. 

in. 

in. 


in. 

Ho 

1 .4 

1 1 0 

570 

5.10 

66.1 

17.8 

226 

0.171 

0.0769 

2.23 

0.0222 

1.00 

1.64 

1 1 5 

573 

4.99 

78.9 

15.3 

174 

0.201 

0.0797 

2.52 

0.01 88 

1.01 

2 

34 

540 

1.39 

92.4 

1 1.5 

126 

0.302 

0.0958 

3.15 

0.01 85 

1.01 

2 

135 

580 

5.01 

92.4 

1 1.5 

126 

0.249 

0.0804 

3.10 

0.0143 

0.99 

3 

56 

538 

1.38 

112.0 

4.43 

69 

0 548 

0.1030 

5.32 

0.0090 

0.99 

3 

150 

579 

3.32 

112.0 

4.43 

69 

0.476 

0.0903 

5.27 

0.0075 

0.97 

4 

90 

558 

1.23 

1 1 8.9 

1.69 

52 

0.883 

0.1050 

8.42 

0.0050 

0.99 

4 

255 

580 

3.29 

1 18.9 

1.69 

52 

0.765 

0.0924 

8.28 

0.0041 

0.97 

5 

330 

605 

2.45 

l 1 8.4 

0.715 

52 

1.136 

0.0926 

12.27 

0.0028 

0.97 

a Porame 

ters used in Eq 

i. £5). 










b From Eq. (4). 











^ Estimated effective average value throughout nozzle; from Fig. 6. 
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Fig. 6. Normalized values of the measured boundary-layer parameter, H, in the JPL 20-in. SWT 


in the test section are too low, ranging from 15% low at 
M — 1.4 to 7% low at M = 5 for the high supply-pressure 
cases. It seems reasonable to expect that the comparisons 
(of the calculated values of boundary-layer thicknesses 
with the measured values) should be essentially the same, 
at each nozzle contour, for all values of supply pressure 
as long as the boundary layer is turbulent from the throat 
on. That is, any shortcomings in the application of the 
boundary -layer, integral momentum equation should 
apply equally as well to the high supply-pressure case as 
to the low supply-pressure case. As this is not the situa- 


tion, and since the calculations for the low supply-pressure 
cases were in considerably better agreement with the 
measured values than they were in the high-prssure cases, 
the use of the experimental data of Fig. 1, which de- 
creases with increasing Reynolds number, does not appear 
to be appropriate for use in wind-tunnel boundary-layer 
calculations. Perhaps the problem with these data is even 
more basic: they should not have been assumed to be 
Reynolds-number dependerit (although, in general, they 
do appear to exhibit such a trend), but rather a best- 
level, straight-line fairing should have been made. 
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It becomes clear that more consistent calculations would 
result if the ratio of C f /CV did not drop off with increas- 
ing Reynolds number. For this reason, the nozzle bound- 
ary-layer calculations were repeated using Eq. (2) [rather 
than Eq. (3)] with Eqs. (5) and (6), the results being 
listed in Table 3. For the high supply-pressure case, the 


Table 2. Calculated values of test-section 
boundary-layer thicknesses using Eq. (3) 



B 

e 

in. 

B 

e 

6’ 

1.4 

High 

0.0653 

0.144 

0.850 

0.843 

1.64 

High 

0.0680 

0.174 

0.853 

0.866 

2 

Low 

0.0917 

0.290 

0.958 

0.962 

2 

High 

0.07)5 

0.222 

0.888 

0.892 

3 

Low 

0.0995 

0.536 

0.966 

0.978 

3 

High 

0.0841 

0.444 

0.932 

0.934 

4 

Low 

0.101 8 

0.873 

0.970 

0.988 

4 

High 

0.0844 

0.709 

0.914 

0.928 

5 

High 

0.0857 

1.065 

0.926 

0.938 



Fig. 7. Analytic value of the turbulent boundary- 
layer parameter, H 0 


Table 3. Calculated values of test-section 
boundary-layer thicknesses using Eq. (2) 


M 

Pt 

e 

in. 

< 5 * 

in. 

e 

5 * 

1.4 

High 

0.071 1 

0.157 

0.924 

0.918 

1.64 

High 

0.0744 

0.191 

0.934 

0.950 

2 

Low 

0.0966 

0.306 

1.008 

1.013 

2 

High 

0.0785 

0.244 

0.976 

0.980 

3 

Low 

0.1051 

0.566 

1.020 

1.033 

3 

High 

0.0914 

0.483 

1.012 

1.015 

4 

Low 

0.1072 

0.919 

1.021 

1.041 

4 

High 

0.0918 

0.771 

0.993 

1.008 

5 

High 

0.0918 

1.140 

0.992 

1.004 


calculated values of the boundary-layer thicknesses range 
from 8% low at M — 1.4 to 2% high at M — 3 in com- 
parison with the measured values. Although the range of 
disagreement from the measured values is not improved 
over that when Eq. (3) was used, at least it now brackets 
unity -a very desirable result. In addition, the effect of 
Reynolds number upon the comparison has been con- 
siderably reduced. Certainly it would be possible to come 
up with a C,/C,. vs. R, relationship that would make the 
calculated values of boundary-layer thickness agree al- 
most exactly with the measured values, but for all practi- 
cal purposes, the- comparisons by use of Eq. (2) are 
adequate. Only at the two lower nozzle Mach numbers 
(M = 1.4, 1.64) are the discrepancies large (5 to 8% ), 
and it is here that they are least important because the 
boundary-layer displacement thickness is such a small 
part of the nozzle height in comparison to the higher 
Mach- number nozzles. 

The Cr/Cq. values of Eq. (2) are compared with those 
of Ref. 12 in Fig. 8. The comparison is quite favorable 
for M < 3, but at M — 6 there is a rather large discrep- 
ancy of about 25%. It is not the intent of this Report to 
discuss the validity of one Reference over another, but 
rather to show that care must be used in choosing the 
local skin-friction law for the turbulent boundary layers 
over an adiabatic flat plate when calculating wind-tunnel 
nozzle boundary layers. 

Figures 9 and 10 indicate the portions of the C f /C f . vs. 
R e regions that are used in the calculation of the 
boundary-layer growth along the nozzle. Upon compari- 
son of these two figures with Fig. 1 it becomes apparent 
that for the high supply-pressure case, the major portion 
of the boundary-layer growth takes place beyond the 
Re region of experimental values of C f /C f . while the 
lower pressure case is in the region of experimental, local 
skin-friction measurements. 
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= 2 


SOLID FAIRINGS FROM EQUATION (2) AND Ref. 10 
OASHED FAIRINGS FROM FlG.7 OF Ref. 12 


. 8. Comparison of the local skin-friction coefficients; in Refs. 10 and 12, of the turbulent boundary layer 

over an adiabatic flat plate 



«• (H- m * J = 1.40 NOZZLE 

m " " = 1.64 NOZZLE 


M -- 3 NOZZLE 


M = 4 NOZZLE 


M -- 5 NOZZLE 
^ ih 


4? = 2 NOZZLE 


0 THROAT 

1 1/4 OF DISTANCE FROM THROAT TO END OF CURVED WALL 

II l/2 OF DISTANCE FROM THROAT TO END OF CURVED WALL 

III 3/4 OF DISTANCE FROM THROAT TO END OF CURVED WALL 

till END OF CURVED WALL 

* POSITION OF BOUNDARY LAYER MEASUREMENTS 


Fig. 9. Cf/Cf. vs. R„ regimes for the calculations of the curved-wall boundary layer for the JPL 20-in. 

SWT nozzle (high supply-pressure cases) 
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Fig. 10. Cf/Cf, vs. R« regimes for the calculations of the curved-wall boundary layer for the JPL 20-in. 

SWT nozzle (low supply-pressure cases) 


The results of the boundary-layer calculations along the parison of these two figures (and Figs. 13 and 14) with Fig. 

nozzle using Eqs. (2), (5), and (6) are shown in Figs. 11 5 shows that the major portion of the boundary-layer 

and 12 and are compared with the measured values in the growth for each nozzle takes place at Mach numbers very 

test section at the high supply-pressure condition. Com- much closer to the test-section Mach number than to the 



X % in. 


Fig. 1 1 . Calculated values of the boundary-layer momentum thickness along the curved nozzle wall of 

the JPL 20-in. SWT (high supply-pressure cases) 
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Fig 12. Calculated values of the boundary-layer displacement thickness along the curved nozzle wall of 

the JPL 20-in. SWT (high supply-pressure cases) 


throat Mach number of M — 1. The effect of the supply 
pressure (high vs. low) upon the boundary-layer growth 
of the M = 3 nozzle is shown in Figs. 13 and 14. 


Figure 15 shows the value of 8 0 [as calculated with 
Eq. (4)] used for both the high and low supply -pressure 
cases. The effect on the test-section boundary-layer mo- 



Fig. 13. Calculated values of the boundary-layer momentum thickness along the M — 3 curved nozzle 
wall of the JPL 20-in. SWT (high vs. low supply-pressure cases) 
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Fig. 14. Calculated values of the boundary-layer displacement thickness along the M =3 curved nozzle 
wall of the JPL 20-in. SWT (high vs. low supply-pressure cases) 


mentum thickness ( at the location of the measured v alues ) 
of the chosen value for the throat momentum thickness 



Fig. 15. Calculated throat boundary-layer momentum- 
thickness as a function of Mach number 


is shown in Fig. 16 for the M - 1.4 and M = 3 nozzles 
(high supply-pressure case). At M - 1.4 a 50% decrease 
or inciease in 0 () results in about a D% decrease or 
increase in the test section 0, while at M = 3 a 100% 
increase in ft, results in only a 1 Vi % increase in the test 
section 0. Figure 1/ presents the proportion change in 
test section 0 relative to a change of one ft, in throat mo- 



THROAT Q 

Fig. 1 6. Effect of chosen value of nozzle throat boundary- 
layer momentum thickness upon calculated value 
of test-section boundary-layer momentum 
thickness (high supply-pressure cases) 


to 
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men turn thickness at <9„. Although theoretical data were 
obtained at only two nozzle Mach numbers, some liberty 
was taken in fairing these two points in order to demon- 


strate the decreasing effect of the chosen 8 U on the test 
section 9 which accompanies an increasing nozzle Mach 
number. 



Fig. 17, Proportional effect of chosen value of 
nozzle-throat boundary-layer momentum 
thickness upon calculated value of test- 
section boundary-layer momentum 
thickness as a function of Mach 
number (high supply- 
pressure cases) 
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IV. CONCLUSIONS 


An analytical formulation was made of experimental, 
local-skin-friction data from the turbulent boundary layer 
over an adiabatic flat plate and was then applied to the 
boundary-layer growth along the curved walls of the JPL 
20-in. Supersonic Wind Tunnel using the boundary-layer 
integral-momentum equation. The resulting test-section 
boundary-layer thicknesses were not consistent with meas- 
ured thicknesses, being generally too small. However, 
the major problem was the variation with Reynolds num- 
ber of the discrepancy between the calculated and 
measured values of the boundary-layer thickness as the 
Reynolds number was changed by a factor of three. 
This was true for several nozzle shapes. 

In order to decrease this variation in discrepancy with 
Reynolds number, the simpler skin-friction formulation 
of Ref. 10 was used in the calculation of the boundary- 
layer growth along the nozzle. In addition, the agreement 
with measured values was considerably improved. For 
2 < M -<5, the calculated values of the boundary-layer 
momentum and displacement-thicknesses agree with the 
measured values to within 4%, even for variations of a 
factor of three in the Reynolds numbers. At the lowest 
Mach number of 1.4, the agreement is still within 8% — 
more than satisfactory for nozzle design. 

The simple, adiabatic, flat-plate, turbulent-boundary- 
layer, local-skin-friction formula of Ref. 10 was integrated 
to give values of R g vs. R x for a flat plate. It was then 


possible to derive an equation of R e vs. R r good to 2% 
for 10 6 < R z < 10 9 and 0 < M < 6 (Appendix A). This 
equation then can be used to estimate the turbulent- 
boundary-layer momentum thickness on an adiabatic flat 
plate. 

A comparison between the flat-plate boundary-layer 
(based upon the R g vs. Zb equation, distance from nozzle 
throat to test section, and test-section unit Reynolds num- 
ber) and the actual measured values in the JPL 20-in. 
SWT can be used to estimate test-section boundary 
layers on the curved walls in two-dimensional wind tun- 
nels of various sizes (see Appendix B). Such a procedure 
gives boundary -layer momentum thicknesses that are 
generally within 5% of the measured values for the three 
tunnels investigated: the JPL 12-in. supersonic wind 
tunnel, with both 9 X 12 and 12 X 12-in. test sections; 
the AEDC 40-in. SWT; and the CIT-CWT 8Mi X HM-ft 
test section. The use of the JPL 20-in. SWT-measured 
ratio of boundary-layer displacement to momentum thick- 
ness gave corresponding values of displacement thick- 
nesses of these three other tunnels which agree to within 
about 13% of the measured values. The tunnel scale 
appears to have a large effect on the 8*/6 ratio even for 
the same value of R x . The data of this investigation indi- 
cated a scale effect, and perhaps these data can be used 
for other tunnels to a higher degree of accuracy than is 
possible by using only the results from the JPL 20-in. 
SWT. 


12 
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APPENDIX A 

R„ vs. R x Relationships 


In many adiabatic, flat-plate, turbulent-boundary-layer 
investigations, the relationship of the boundary-layer 
momentum-thickness Reynolds number ( ZL ) as a function 
of the flat-plate-length Reynolds number ( R. r ) is required. 
This relation can be obtained by integrating Eq. (A-l) 
using the experimental data skin friction formula of Eq. 
(3) for M = 0, 1, 2, 3, 4, 5, and 6. 



The distance increments (Ax) used for the integration of 
Eq. (A-l) were: 


x, in. 

Ax, in. 

0 to 5 

0.1 

5 to 10 

0.5 

10 to 50 

1 

50 to 100 

5 

100 to 500 

10 

500 to 1000 

50 

1000 to 5000 

100 

5000 to 10000 

500 


With a unit Reynolds number of 10" per in., the value of 
the momentum thickness (0,,) at the plate leading edge 
(x — 0) was taken to be 0.001 in. The results of this inte- 
gration are summarized in Table A-l. 

The effect of varying the leading-edge momentum 
thickness (0„) was investigated by letting = 0.003 in., 
and is shown in Table A-2. At R r — 10 (> , the variation in 
R 0 due to a change in 0„ from 0.001 in. to 0.003 in. is 
about 5%, and at ZT — 10 7 it is down to a negligible 
Vi%. Therefore, the choice of is not very important 
as long as it is small (around 0,001 in,). The effect of 
varying the size of the Ax increment lengths was investi- 
gated by the expedient procedure of decreasing the unit 
Reynolds number from 10 to 5 X 10 4 per in. This, in 
essence, halved the Ax increment lengths at a unit Rey- 
nolds number of 10 . The resulting values of R 0 as a 
function of R, (summarized in Table A-2) indicate an 
effect of only 1% at fl, — 10' ; due to smaller incremental 
lengths along the flat plate in the integration of Eq. 
(A-l). 


The same integration process used in Eq. (A-l) was 
performed using the skin-friction formula of Eq. (2) 
(see Fig. A-l ) which is independent of Reynolds number. 
(Although the ratio of C//C/. can be assumed independ- 
ent of Reynolds number, C/. is strongly dependent upon 
Reynolds number.) The results of this integration are 
summarized in Table A-3 and in Fig. A-2. 

As the adiabatic flat-plate, turbulent-boundary-layer 
relationship of Fig. A-2 ( R 0 vs. R, ) is quite useful, an 
attempt was made to obtain an analytical formulation of 
ZE as a function of R x in order that interpolation of the 
data could be simplified while retaining a fair degree of 
accuracy. Based upon the skin-friction law of Eq. (2), 
the analytical expression relating R» to R r (based upon 
the integration procedure discussed previously in this 
Appendix) is: 

R e = A X B X C (A-2) 

where 

A is the M = 0 term 

B is the M effect term at R x = 10“ 

C is the Rj effect on the R 0 term for M > 0 and R x ^ 1(1’ 



Fig. A-l. Cf/C/. vs. M 
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The actual algebraic forms for A, B, and C are: 


A = 

0.226 fl. 

(A-3) 

[i 

1 0.033 (Ioe w R, - 6.7)] [log U) R,]*-** 

B ( 1 

f 0.144 M 2 ) n 4,1 

(A-4) 


M 

l- 1 

t —(Ioz,„R,. - 6) (1 -1- 0.144 M 2 )<'•« 

"J 

(*-5) 


The degree to which Eq. ( A-2) fits the data of Table A -3 
and Fig. A-2 is indicated in Table A -4 where the devia- 
tion of the Eq. (A-2) values from the integrated values 
of Ro vs. R r is indicated in per cent for 0 < M < 6 and 
5 X 10 :> < Rr < 10”. For R; ^ 10' the errors are gener- 
ally less than 2%, a more 1 than satisfactory match of the 
data by an equation. 


The “A” term of Eq. (A-2) is merely an alteration of 
the Schlichting M - 0 relation [Ref. 12; see Eq. (A-6)] 
in order to obtain a better fit of the integrated values of 
the Schultz-Crunow relation of Eq. (1 ). 


C f j 


0.455 

(log i„ R ,) 2 r,s 


(A-6) 



Fig. A-2. Ro vs. R x as a function of Mach number 
[see Eq. (A-2)] 


For M 0 only, the constant term of Eq. (A-3) should the constant of Eq. (A-6). The “B” term of Eq. (A-2) is 

be increased by about J% yielding one-half the value of simply a minor modification of Eq. (2). 
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Table A-l . Adiabatic flat-plate, turbulent-boundary-layer properties 
using Eq. (3); R/in. = 10 s , (9 0 = 0.001 in. 
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Table A-3. Adiabatic flat-plate turbulent-boundary-layer properties 
using Eq. (2); R/in.= 10 s , = ' n - 
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Table A-4. Percent error" of R 0 vs. R x of Eq. (A-2) in matching 
the data of Table A-3 


** 

M — 0 

M = 1 

M — 2 

M = 3 

M = 4 

M = 5 

M = 6 

5 X 10 5 

- 4.5 

- 4.0 

- 4.5 

- 4.4 

- 5.2 

- 6.2 

- 6.2 

1 X 10 4 

- 1.7 

- 1.7 

- 1.8 

- 1.6 

- 1.6 

- 1.7 

- 1.8 

2 X 10 * 

- 0.9 

- 0.5 

- 0.4 

- 0.2 

- 0.1 

0 

- 0.2 

5 X 10 * 

0.2 

0.1 

0.8 

1.3 

1.4 

1.5 

1.1 

1 X 10 7 

- 0.7 

- 0.2 

0.8 

0.7 

1.2 

1.2 

1.4 

2 X 10 7 

0.7 

0.3 

1.1 

1.9 

1.8 

2.3 

1.8 

5 X 10 7 

1.0 

0.2 

0.4 

1.3 

1.6 

1 .1 

0.4 

1 X 10 8 

- 0.5 

- 0.1 

0.5 

0.9 

1.3 

1.2 

- 0.2 

2 X 10 * 

- 0.9 

- 0.9 

- 0.9 

- 0.1 

0.1 

- 0.3 

- 1.2 

5 X 10 8 

- 1.7 

- 1.6 

- 0.9 

- 0.6 

- 0.9 

- 1.3 

- 2.6 

1 X 10 ? 

- 1.3 

- 1.6 

- 1.1 

- 0.7 

- 0.9 

- 1.4 

- 3.0 

“These percent error values are approximate, having been calculated to slide rule accuracy only. 


APPENDIX B 

Short-Cut Method of Estimating Test-Section Boundary-Layer 
Momentum and Displacement Thickness 


The comparisons of the calculations of the test-section 
boundary layer along the curved wall of the JPL 20-in. 
SWT with the turbulent boundary-layer growth along 
an adiabatic flat plate can be used to estimate the corre- 
sponding boundary layers of other two-dimensional wind 
tunnels. Table B-l compares the adiabatic flat-plate, tur- 
bulent boundary layer (based upon the distance from 
the throat to the boundary-layer station in the test sec- 
tion, and assuming test-section Mach number and 
Reynolds number throughout) to the measured values in 
the test section. These ratios are plotted in Fig. B-l 
along with similar results from several other facilities, 
both smaller and larger. 


The ratios of Fig. B-l in conjunction with the H/H 0 
values of Fig. 6 can be combined to give rough estima- 
tions of the curved-wall, turbulent-boundary-layer mo- 
mentum and displacement thickness in the test section. 
The flat-plate boundary-layer data of Table A-3, Fig. A-2 
or Eq. (A-2) are used for the ratios of Fig. B-l. The 
results of these calculations are shown in Tables B-2, B-3, 
and B-4 for the JPL 12-in. SWT (Ref. 8), CIT Co-op 
Wind Tunnel (Ref. 13), and the AEDC 40-in. SWT 
( Refs. 14, 15 ) . These boundary-layer momentum-thickness 
estimations are generally good to within about 5%. 
However, the displacement thicknesses based upon this 
simple method are off by as much as 13%, and, in gen- 
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era), are never better than about 5%. The reasons for 
these discrepancies can be seen from Figs. B-l and B-2. 

The relationship between di,, and 9 measurements in 
Fig. B-l is not monotonic with tunnel size, as the JPL 
20-in. SWT curve is straddled by the curves obtained 
from larger facilities. Comparing the data as a function 



Fig. B-l. Ofp/Om vs. M for several supersonic 
wind tunnels 


of tunnel test-section R , does not offer any more hope 
than the use of tunnel size. Perhaps a universal (average) 
curve would be the most suitable. It is quite definite that 
in a given facility, the curve moves upward with decreas- 
ing Reynolds number. 

The ratios of H/H„ in Fig. B-2 appear to be monotonic 
with tunnel size, and in a given facility, the curve is 
raised with a decrease in Reynolds number. Perhaps a 
better estimate of 8* can be made by allowing for the 
effect of tunnel size than was obtained by simply using 
the results from the JPL 20-in. SWT shown in Fig. 6. 

The basic relationships of the short-cut method are as 
follows: Ro is a function of R, and M [Eq. (A-2) or Fig. 
A-2]. Use M for M. R x is based on: (a) distance from 
nozzle throat to test-section Jpoundary-layer station; (b) 
test-section Mach number (A/); and (c) test-section unit 
Reynolds number (R/in.). 

0,p R/in. 


Oailr — fyp j q 

§*calc = 9 C ah/ X H„ X JJ- 

The value of 9 u ,/9 m is from the JPL 20-in. SWT curve 
in Fig. B-l. The value of HJlh is from the JPL 20-in. 
SWT curve in Fig. B-2. 



Fig. B-2. Normalized values of the measured boundary-layer parameter, H, from several supersonic 

wind tunnels 
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Table B-l. JPL 20-in. SWT test-section boundary- layer measurements 
on centerline of curved wall 


M 

m 

M 

#x 

X HT 7 

5*m 

in. 

in. 


*9 X 10 4 

(flat plate) 

in. 

hp. 

e m 

1.40 

5.10 

66,1 

3.37 

0.171 

0.0769 

2.23 

3.59 

0.0705 


1 .64 

4.99 

78.9 

3.94 

0.201 

0.0797 

2.52 

3.94 

0.0785 

■ 

2.01 

1.39 

92.4 

1.28 

0.302 

0.0958 

3.15 

1.46 

0.1050 


2.01 

5.01 

92.4 

4.63 

0.249 

0.0804 

3.10 

4.24 

0.0848 


2.61 

3.65 

105.2 

3.84 

0.371 

0.0868 

4.27 

3.23 

0.0885 

1.021 

2.61 

3.70 


3.84 

0.364 

0.0845 

4.31 

3.23 

0.0873 

1.033 

3.00 

1.38 

11 2.0 

1.54 

0.548 

0.1030 

5.32 

1.42 

0.1029 

0.998 

3.00 

3.32 

1 12.0 

3.72 

0.476 

0.0903 

5.27 

2.95 

0.0889 

0.985 

3.50 

3.34 

1 17.1 

3.91 

0.603 

0.0899 

6.71 

2.80 

0.0838 

0.933 

3.74 

3.1 8 

1 18.3 

3.76 

0.660 

0.0866 

7.63 

2.59 

0.0815 

0.942 

3.98 

1.23 

1 18.9 

1.46 

0.883 

0.1050 

8.42 

1.14 

0.0927 

0.882 

3.98 

3.29 

1 18.9 

2.90 

0.765 

0.0924 

8.28 

2.58 

0.0784 

0.848 

4.23 

3.54 

11 8.9 

4.21 

0.823 

0.0891 

9.23 

2.65 

0.0749 

0.841 

4.54 

3.36 

1 18.9 

4.00 

0.983 

0.0944 

10.41 

2.41 

0.0717 

0.760 

4.76 

3.02 

1 18.6 

3.58 

1.052 

0.0921 

1 1.42 

2.1 2 

0.0703 

0.764 

5.00 

2.45 

1 18.4 

2.90 

1.136 

0.0926 

12.27 

1.72 

0.0703 

0.760 


Comparison ^ 

« 

with flat-plate boundary layer. 


Table B-2. JPL 12-in. SWT test-section boundary-layer measurements 
on centerline of curved wall 





R x 



H m 

* e X 10-< 



&\alc 

e calc 

5 *calc 



in. 

X 10~ 7 

in. 

in. 


(flat plate) 


1 

in. 

e m 


|E| 

5.62 

36.5 

2.05 

BH 

0.0517 

2.07 

2.42 

0.0431 

0.0487 

0.099 

0.94 

0.93 

mm 

5.77 

42.5 

2.45 

1 19 

0.0539 

2.28 

2.76 

0.0478 

0.0510 

0.1 15 

0.95 

0.94 

1.61 

5.92 

47.5 

2.81 

■SB 

0.051 2 

2.62 

2.98 

0.0504 

0.0515 

0.127 

1.01 

0.95 

1.78 

5.93 

50.5 

2.99 

0.150 

0.0518 

2.90 

3.06 

0.0516 

0.0504 

0.138 

0.97 

0.92 

1.97 

5.51 

53.0 

2.92 

0.165 

0.0520 

3.18 

2.90 

0.0527 

0.0507 

0.152 

0.98 

0.92 

2.1 8 

5.29 

55.0 

2.91 

0.186 

0.0515 

3.62 

2.77 

0.0524 

0.0504 

0.172 

0.98 

0.92 

2.20 

6.20 

45.9 

2.85 

0.153 

0.0423 

3.62 

2.73 

0.0440 

0.0423 

0.145 

1.00 

0.95 

2.36 

6.38 

47.0 

3.00 

0.170 

0.0427 

3.98 

2.77 

0.0434 

0.0417 

0.157 

0.98 

0.92 

2.54 

6.54 

47.9 

3.13 

0.184 

0.0415 

4.44 

2.75 

0.0421 

0.0408 

0.169 

0.98 

0.92 

2.70 

6.7 8 

48.7 

3.30 

0.205 

0.0424 

4.84 

2.82 

0.0416 

0.0407 

0.182 

0.96 

0.89 

2.80 

6.75 

49.1 

3.31 

0.213 

0.0419 

5.09 

2.77 

0.041 1 

0.0403 

0.191 

0.96 

0.90 

2.98 

6.87 

49.7 

3.41 

0.229 

0.0435 

5.27 

2.77 

0.0403 

0.0409 

0.213 

0.94 

0.93 








Estimation of test-section boun 

idary layer u 

sing results 

of 









JPL 20 in.-SWT measurements. 
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Table B-3. CWT test-section boundary-layer measurements on 
centerline of curved wall 


M 


mm 

*x 


0m 

H m 

R«X 10- 1 


0c ale 

5%a/c 

^ cot c 

5 \a/e 



X icr 7 

in. 

in. 


(flat plate) 

in. 

in. 

in. 

0m 

5 m 

1.21 

■ 


3.86 

0.519 

, 0.284 

1.83 

4.14 

0.258 

0.293 

0.588 

1.03 

1.13 

1.30 

1.40 

240 

3.37 

0.515 

0.261 

1.98 

3.63 

0.259 

0.270 

0.558 

1.03 

1.08 

1.40 

1.49 

240 

3.58 

0.531 

0.252 

2.1 1 

3.77 

0.253 

0.263 

0.597 

1.04 

1.12 

1.50 

1.48 

240 

3.55 

0.541 

0.243 

2.22 

3.68 

0.248 

0.248 

0.595 

1.02 

1.10 

1.60 

1.32 

240 

3.1 8 

0.548 

0.232 

2.36 

3.30 

0.250 

0.245 

0.617 

1.06 

1.13 

1.70 

1.33 

240 

3.20 

0.558 

0.222 

2.51 

3.27 

0.246 

0.239 

0.633 

1.08 

1.13 

1.75 

1.29 

240 

3.10 

0.572 

0.224 

2.55 

3.16 

0.245 

0.236 

0.647 

1.05 

1.13 






Estimation of test-section boundary layer using results of 
JPL 20-in. SWT measurements. 


Table B-4. AEDC 40-in. SWT test-section boundary-layer measurements 
on centerline of curved wall 


a 



X 10 7 

£ m 

in. 

0m 

in. 

" n , 

X 10* 

(flat plate) 

% 

in. 

0ca/c 

in. 


0co/c 

0m 

<5*co/c 


1.96 

214.5 

4.21 

0.415 

0.1 80 

2.31 

4.25 

0.181 

0.189 

0.442 

1.05 

1.06 

2 

3.50 

214 

7.48 

0.475 

0.160 

2.97 

6.37 

0.152 

0.148 

0.467 

0.93 

0.98 

3 

1.75 

215 

3.76 

0.92 

0.180 

5.1 1 

29.7 

0.170 

0.170 

0.915 

0.94 

1.00 

3 

5.25 

21 5 

1 1.30 

0.78 

0.155 

5.03 

7.51 

0.143 

0.144 

0.761 

0.93 

0.98 

4 

0.93 

215 

2.01 

1.54 

0.190 

8.1 1 

1.49 

0.159 

0.177 

1.52 

0.93 

0.99 

4 

5.00 

215 

10.70 

1.1 8 

0.148 

7.97 

6.00 

0.120 

0.139 

1.17 

0.94 

0.99 

5 

5.67 

214 

12.10 

1.65 

i 

0.137 

12.04 

5.67 

0.100 

0.135 

1.66 

0.99 

1.00 



These meal 
re! 

turements ej 
suits in Refs 

itimated from plotted 
14 and 15. 

Estimation of test-section boundary layer using results of 
JPL 20-in. SWT calculations. 
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NOMENCLATURE 


Cj compressible local skin-friction coefficient 
Cf incompressible local skin-friction coefficient 

h a */0 

H„ nominal H = 1.3 + 0.46M 2 

L* nozzle throat radius of curvature 

M Mach number: 1) along nozzle 

2) over flat plate 

M test-section Mach number 

P, supply-section pressure 

R Reynolds number 

R/in. Reynolds number per inch 

R x length Reynolds number : xR/ in. 

R„ momentum Reynolds number — OR/ in. 

SWT supersonic wind tunnel 

t* nozzle throat height 

T, supply-section temperature 


x distance from: 1) throat of nozzle 

2 ) leading edge of flat plate 

Ax incremental value of x 

8* boundary-layer displacement thickness 

6 boundary-layer momentum thickness 

0 fp turbulent boundary-layer momentum thickness 
for adiabatic flat plate 

0 a boundary-layer momentum thickness at nozzle 
throat 

e ts test-section boundary-layer momentum thick- 
ness as calculated by boundaiy -layer 
momentum equation 

v* air viscosity in nozzle-throat boundary layer 
based on average of wall and freestream 
temperatures 

( ), alc calculated values using short-cut method of 

Appendix B 

( ) m measured values 
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